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^ Abstract 
CO 

We derive compact analytical formulae for all tree-level color-ordered gauge the- 
ory amplitudes involving any number of external gluons and up to four massless 
quark-anti-quark pairs. A general formula is presented based on the combinatorics 
> of paths along a rooted tree and associated determinants. Explicit expressions 

are displayed for the next-to-maximally helicity violating (NMHV) and next-to- 
next-to-maximally helicity violating (NNMHV) gauge theory amplitudes. Our 
results are obtained by projecting the previously-found expressions for the super- 
amplitudes of the maximally supersymmetric super Yang-Mills theory (Af = 4 
SYM) onto the relevant components yielding all gluon-gluino tree amplitudes in 
M = 4 SYM. We show how these results carry over to the corresponding QCD 
amplitudes, including massless quarks of different flavors as well as a single elec- 
troweak vector boson. The public Mathematica package GGT is described, which 
encodes the results of this work and yields analytical formulae for all Af = 4 SYM 
gluon-gluino trees. These in turn yield all QCD trees with up to four external 
arbitrary-flavored massless quark-anti-quark pairs. 
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1 Introduction and conclusions 

Scattering amplitudes play a central role in gauge theory. At a phenomenological level, they 
are critical to the prediction of cross sections at high-energy colliders, for processes within and 
beyond the Standard Model. Efficient evaluation of scattering amplitudes involving many quarks 
and gluons is particularly important at machines such as the Large Hadron Collider (LHC), in 
which complex, multi-jet final states are produced copiously and complicate the search for new 
physics. 

Tree amplitudes can be used to predict cross sections at leading order (LO) in the perturbative 
expansion in the QCD coupling a s . Such results are already available numerically for a wide 
variety of processes. Programs such as MadGraph |1|, CompHEP |2|, and AMEGIC + + |3| 
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are based on fast numerical evaluation of Feynman diagrams. Other methods include the Berends- 
Giele (off shell) recursion relations |4|, as implemented for example in COMIX |5|, and the 
related ALPHA |6| and HELAC |7| algorithms based on Dyson- Schwinger equations, as well as 
O' Mega/ WHIZ ARD |8|. The computation time required in these latter methods scales quite 
well with the number of legs. 

On the formal side, the properties of scattering amplitudes have long provided numerous 
clues to hidden symmetries and dynamical structures in gauge theory. It was recognized early 
on that tree amplitudes in gauge theory are effectively supersymmetric |9|, so that they obey 
supersymmetric S- matrix Ward identities |10|. Soon thereafter, Parke and Taylor |11| discovered 
a remarkably simple formula for the maximally-helicity-violating (MHV) amplitudes for n-gluon 
scattering, which was proven by Berends and Giele |4|, and soon generalized to M = 4 super 
Yang-Mills theory (SYM) by Nair [12]. 

Later, it was found that this simplicity also extends to the loop level, at least for M = 4 
super Yang-Mills theory |13, 14 1. These results were obtained using the unitarity method, which 
constructs loop amplitudes by sewing together tree amplitudes (for recent reviews see refs. |15|). 
After Witten |16| reformulated gauge theory in terms of a topological string propagating in 
twistor space, there was a huge resurgence of interest in uncovering new properties of scattering 
amplitudes and developing new methods for their efficient computation. Among other devel- 
opments, Britto, Cachazo, Feng and Witten proved a new type of recursion relation |17| for 
gauge theory. In contrast to the earlier off-shell recursion relations, the BCFW relation uses only 
on-shell lower-point amplitudes, evaluated at complex momenta. A particular solution to this 
recursion relation was found for an arbitary number of gluons in the split-helicity configuration 
( +■••+) [18]. 

The BCFW recursion relation was then recast as a super-recursion relation for the tree am- 
plitudes of M = 4 super Yang-Mills theory, which involves shifts of Grassmann parameters as 
well as momenta |19|. A related construction is given in ref. |20|. The super-recursion rela- 
tion of ref. 1 19 1 was solved for arbitrary external states by Drummond and one of the present 
authors |21|. Tree-level super-amplitudes have a dual superconformal invariance |22,23|, and 
the explicit solution does indeed have this symmetry |21|. It is written in terms of dual super- 
conformal invariants, which are a straightforward generalization of those that first appeared in 
next-to-MHV (NMHV) super-amplitudes |22,24|. This dual superconformal invariance of tree- 
level amplitudes is a hallmark of the integrability of planar = 4 SYM, as it closes with the 
standard superconformal symmetry into an infinite-dimensional symmetry of Yangian type |25| 
(a recent review is ref. |26|). 

The purpose of this paper is to illustrate how these more recent formal developments can reap 
benefits for phenomenological applications in QCD. In particular, we will evaluate the solution 
in ref. |21| by carrying out the integrations over Grassmann parameters that are needed to select 
particular external states. In addition, we will show how to extract tree-level QCD amplitudes 
from the amplitudes of Af = 4 super Yang-Mills theory. While this extraction is simple for pure- 
gluon amplitudes, and those with a single massless quark line, it becomes a bit more intricate 
for amplitudes with multiple quark lines of different flavors, because of the need to forbid the 
exchange of scalar particles, which are present in M = 4 super Yang-Mills theory but not in 
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QCD. 

Although, as mentioned above, there are currently many numerical programs available for 
computing tree amplitudes efficiently, the existence of analytic expressions may provide a yet 
more efficient approach in some contexts. In fact, the formulae provided in this paper have 
already served a practical purpose: They were used to evaluate contributions from real emission 
in the NLO corrections to the cross section for producing a W boson in association with four 
jets at the LHC |27|. This process forms an important background to searches for various kinds 
of new physics, including supersymmetry. The real-emission corrections require evaluating nine- 
point tree amplitudes at a large number of different phase-space points (on the order of 10 8 ), in 
order to get good statistical accuracy for the Monte Carlo integration over phase space. 

In principle, QCD tree amplitudes can also be used to speed up the evaluation of one-loop 
amplitudes, when the latter are constructed from tree amplitudes using a numerical implementa- 
tion of generalized unitarity. Many different generalized unitarity cuts, and hence many different 
tree amplitudes, are involved in the construction of a single one-loop amplitude. The tree ampli- 
tudes described here enter directly into the construction of the "cut-constructible" part |14| of 
one- loop amplitudes in current programs such as CutTools |28|, Rocket |29,30| and Black- 
Hat 1 3 X | . On the other hand, the computation-time bottleneck in these programs often comes 
from the so-called "rational" terms. When these terms are computed using only unitarity, it is 
via unitarity in D dimensions 1 32 , 30 , 33 1 , not four dimensions. The amplitudes presented here 
are four- dimensional ones, so they cannot be used directly to alleviate this bottleneck for the D- 
dimensional unitarity method. However, in the numerical implementation of loop-level on-shell 
recursion relations |34| for the rational part in BlackHat |31|, or in the OPP method used in 
CutTools |28|, there are no .D-dimensional trees, so this is not an issue. 

An interesting avenue for future research would be to try to generalize the results presented 
in this paper to QCD amplitudes containing massive quarks, or other massive colored states. 
Massive quark amplitudes are of interest because, for example, processes that produce top quarks 
in association with additional jets can form important backgrounds to new physics at the LHC. 
States in M = 4 SYM can be given masses through a super-Higgs mechanism. This mechanism 
was explored recently in the context of infrared regulation of M = 4 SYM loop amplitudes |35|. 
However, it should be possible to generate the appropriate tree amplitudes with massive quarks, 
or other massive states, from the same kind of setup, once one solves the appropriate super- 
BCFW recursion relations. 

The remainder of this paper is organized as follows. After introducing the standard technology 
of color-ordered amplitudes and spinor helicity we explain the strategies of how to extract QCD 
tree amplitudes with massless quarks from M = 4 SYM in section 3. We also discuss how to 
convert these amplitudes into trees with one electroweak vector boson. Sections 4 through 6 are 
devoted to stating the general analytical formulae for gluon-gluino n-parton amplitudes in = 4 
SYM, which are proven in section 7. In appendix A we provide a collection of explicit results for 
pure-gluon trees. Explicit formulae for trees involving up to six fermions are displayed in appendix 
B. Finally, appendix C is devoted to a documentation of our Mathematica package GGT which 
implements all of the results of this paper and yields the analytical expressions for an arbitrary 
flavored gluon-gluino tree amplitude in M = 4 SYM. The package is included in the arXiv.org 
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submission of this article and may also be downloaded from http : //qf t . physik.hu-berlin. de. 



2 Color-ordering and spinor-helicity formalism 

Tree-level gluon amplitudes in non-abelian gauge theories may be conveniently separated into a 
sum of terms, each composed of a simple prefactor containing the color indices, multiplied by a 
kinematical factor known as a partial or color-ordered amplitude. For an n-gluon amplitude one 
has 

A t r r({Pi,h i ,a l })=g n - 2 Tr(nw...r^))4(a(l)^...aH^), (2.1) 

with the argument i h ' of the partial amplitude A n denoting an outgoing gluon of light-like mo- 
mentum pi and helicity hi = ±1, i G [1, n]. The su(N c ) generator matrices T ai are in the 
fundamental representation, and are normalized so that Tr(T a T b ) = 5 ab . 

Color-ordered amplitudes of massless particles are most compactly expressed in the spinor- 
helicity formalism. Here all four-momenta are written as bi-spinors via 



laa 



<V, (2-2) 



where we take = (1, a) with a being the 2x2 Pauli spin matrices. Light-like vectors are then 
expressed via the product of two spinors 

f* = \ a \ & . (2.3) 

For real momenta with Lorentz signature we have A = ±A*, with the sign being determined 
by the energy component of p. For complex momenta the spinors A and A are independent. 
Our convention is such that all gluons are outgoing. Then in eq. (2.1) each color-ordered leg is 
specified by a choice of spinors Aj and Aj along with a helicity hi = ±1. Given this data the 
associated polarization vectors may be reconstructed from the expressions 

\d ,.a \a ~a 

act i r^i aa i r'i (n a\ 

+,i ~ /\ \ ' —>i ~ r\ i ' 
(Ai/ij) [AiHil 

where are auxiliary momenta and we use the standard notation (A/i) = e a p\ a \$ and 

[A//] = e^A"/^. Moreover we shall often use the abbreviated forms (ij) = (AjAj) and [ij] = 
in the sequel. As an essential building block of the general tree-level scattering formula we 
introduce the dual coordinates or region momenta xff via 



j'-i 



to + j> i+x + ■■■+ ^-ir = E xa ^ > * < i> ( 2 - 5 ) 

k=i 

Xu = 0, and Xij = —Xji for i > j. We then define the scalar quantities 

(naia 2 . . . a* | a) := (n\x nai x aia2 . . . x afc _ ia Ja) , (2.6) 
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which we will use frequently in the following. In fact all amplitudes can be expressed in terms of 
the quantities {na\a 2 . . . Ofc|a) and the spinor products (i j). 

As an example of the notation and in order to give a flavor of the kinds of results we obtain, 
we present a compact formula for the n-point NMHV pure gluon amplitude in QCD 



A 



NMHV, 



ii,i2,n) 



6^(p) 



(1 2) ... (n 1) 



^2 R n] st((nii)(rits\i2)^ + ^ n )( n ^i) 37 ^) 



ii<s<i2<t<n— 1 



i\ <s<t< «2 



> i^-atf (^2 ii){nts\n)J + 2^ i2 n;st ( (n z 2 )(nst|zi) 

2<s<ii<J2<i<n-l 2<s<ii<t<i2 



(2.7) 



Here ii,Z2 and n correspond to the positions of the three negative-helicity gluons. (Using cyclic 
symmetry, we have put one of them at position n without loss of generality.) The quantities R n - S t 
are simply given by 



R 



(8(8-1)) 



(t(t-l)) 



n;st 



1 

x 2 st (nts\s) (nts\s — 1) (nst\t) (nst\t — 1) 



(21 



with R n]S t '■= for t = s + lois = t+l. Note that the above formula is given for an arbitrary 
number of gluons n. In realistic cases this number is usually small, say of the order of 9, in which 
case relatively few terms are produced by the nested sums in eq. (2.7). 



3 From AT = 4 SYM to QCD tree amplitudes 

In this section we discuss how to assign quantum numbers for external states in M = 4 SYM in 
order to generate tree amplitudes for QCD with massless quarks. We then discuss the genera- 
tion of tree amplitudes including an electroweak vector boson (W, Z or virtual photon). From 
the point of view of tree amplitudes, there are two principal differences between = 4 SYM 
and massless QCD. First of all, the fermions in = 4 SYM, the gluinos, are in the adjoint 
representation of su(N c ), rather than the fundamental representation, and come in four flavors. 
Secondly, the M = 4 SYM theory contains six massless scalars in the adjoint representation. 

Because we use color-ordered amplitudes, as discussed in Sec. 2, the first difference is fairly 
unimportant. Quark amplitudes can be assembled from the same color-ordered amplitudes as 
gluino amplitudes, weighted with different color factors. For example, the color decomposition 
for amplitudes with a single quark- ant i- quark pair, and the remaining (n — 2) partons gluons is, 

A^ e (l q ,2 q ,3,...,n) = g n - 2 £ (T^ . . ,TM")).J A« ee (h, 2 q , a(3), . . . , a(n)) . (3.1) 

<TS5 n -2 

The color-ordered amplitudes appearing in eq. (3.1) are just the subset of two-gluino-(n— 2)-gluon 
amplitudes in which the two gluinos are adjacent. 
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Amplitudes with more quark-anti-quark pairs have a somewhat more intricate color structure 
involving multiple strings of T a matrices, as explained in ref. |36|. In this case, some of the color 
factors also include explicit factors of l/N c , as required to project out the su(N c ) -singlet state for 
gluon exchange between two different quark lines. However, all of the required kinematical coef- 
ficients can still be constructed from suitable linear combinations of the color-ordered amplitudes 
for 2k external gluinos and (n — 2k) external gluons. 

The main goal of this section will be to illustrate how to choose the flavors of the external 
gluinos in order to accomplish two things: (1) avoid the internal exchange of massless scalars, 
and (2) allow all fermion lines present to be for distinct flavors. (In some cases one may want 
amplitudes with (partially) identical fermions; these can always be constructed from the distinct- 
flavor case by summing over the relevant exchange-terms, although it may be more efficient 
to compute the identical-fermion case directly.) We will accomplish this goal for amplitudes 
containing up to four separate fermion lines, that is, eight external fermion states. 

In gauge theory, tree amplitudes that contain only external gluons are independent of the 
matter states in the theory |9|; hence they are identical between Af = 4 SYM and QCD. The 
reason is simply that the vertices that couple gluons to the other states in the theory always 
produce the fermions and scalars in pairs. There are no vertices that can destroy all the fermions 
and scalars, once they have been produced. If a fermion or scalar were produced at any point 
in a tree diagram, it would have to emerge from the diagram, which would no longer have only 
external gluons. In other words, the pure-gluon theory forms a closed subsector of Af = 4 SYM. 

Another closed subsector of Af = 4 SYM is Af = 1 SYM, which contains a gluon and a single 
gluino. Let g denote the gluon, g^, A = 1,2,3,4, denote the four gluinos, and <Pab = —<Pba 
denote the six real scalars of Af = 4 SYM. Then the Af = 1 SYM subsector is formed by (g, g±). 
The reason it is closed is similar to the pure-gluon case just discussed: There are vertices that 
produce states other than (g,gi), but they always do so in pairs. For example, the Yukawa 
coupling 4> AB gA9Bi A ^ B, can convert g\ into a scalar and a gluino each carrying an index 
B^l. However, this index cannot be destroyed by further interactions. 

The fact that Af = 1 SYM forms a closed subsector of Af = 4 SYM, in addition to color 
ordering, immediately implies that any color-ordered QCD tree amplitude for gluons, plus arbi- 
trarily many quarks of a single flavor, is given directly by the corresponding amplitude (with g\ 
replacing the single quark flavor) evaluated in Af = 4 SYM. The less trivial QCD amplitudes to 
extract are those for multiple fermion flavors, primarily because of the potential for intermediate 
scalar exchange induced by the Yukawa coupling (p AB gAgB- Figure 1 illustrates scalar exchange 
in an amplitude with four fermions belonging to two different flavor lines, A ^ B. 

The key to avoiding such unwanted scalar exchange is provided by figure 2 , which shows four 
types of vertices that could potentially couple fermion pairs to scalars and gluons. However, all 
four types of vertices vanish. (Recall that helicities are labeled in an all-outgoing convention.) 
Case (a) vanishes because the Yukawa interaction only couples gluinos of different flavors, A ^ B. 
Cases (b) and (c) vanish because of fermion helicity conservation for the gauge interactions, and 
a helicity flip for the Yukawa coupling. Case (d) vanishes because gluon interactions do not 
change flavor. 

Because the emission of gluons from fermions does not change their helicity or flavor, in 



7 




Figure 1: Unwanted scalar exchange between fermions of different flavors, A ^ B. 




(a) (b) (c) (d) 

Figure 2: These vertices all vanish, as explained in the text. This fact allows us to 
avoid scalar exchange and control the flow of fermion flavor. 

analyzing whether scalar exchange can be avoided, as well as the pattern of fermion flavor flow, 
one can ignore the giuons altogether. For example, figure 3 shows the possible cases for amplitudes 
with one or two fermion lines. The left-hand side of the equality shows the desired (color-ordered) 
fermion-line flow and helicity assignment for a QCD tree amplitude. All giuons have been omitted, 
and all fermion lines on the left-hand side are assumed to have distinct flavors. The right-hand 
side of the equality displays a choice of gluino flavor that leads to the desired amplitude. All other 
one- and two-fermion-line cases are related to the ones shown by parity or cyclic or reflection 
symmetries. 

The one-fermion line, case (1), is trivial because M = 1 SYM forms a closed subsector of 
M = 4 SYM. In case (2a) we must choose all gluinos to have the same flavor; otherwise a scalar 
would be exchanged in the horizontal direction. Here, helicity conservation prevents the exchange 
of an unwanted gluon in this direction, keeping the two flavors distinct as desired. In case (2b), 
we must use two different gluino flavors, as shown; otherwise helicity conservation would allow 
gluon exchange in the wrong channel, corresponding to identical rather than distinct quarks. 

More generally, in order to avoid scalar exchange, if two color- adjacent gluinos have the same 
helicity, then we should choose them to have the same flavor. In other words, we should forbid all 
configurations of the form (. . . , A + , B + , . . .) and (. . . , A~ , B~ , . . .) for A ^ B, where A ± stands 
for the gluino state g\. While this is necessary, it is not sufficient. For example, we also need 
to forbid configurations such as (. . . , A + , C^, C T , B + , . . .), because the pair (C^, C T ) could be 
produced by a gluon splitting into this pair, which also connects to the (A + ,B + ) fermion line. 
As a secondary consideration, if two color-adjacent gluinos have opposite helicity, then we should 
choose them to have the same flavor or different flavor according to the desired quark flavor flow 
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Figure 3: The possible fermion-line configurations for amplitudes with one fermion 
line, case (1), or two fermion lines, cases (2a) and (2b). 




Figure 4: The possible fermion-line configurations for amplitudes with three fermion 
lines. 



on the left-hand side. However, it will not always be possible to choose them to have a different 
flavor. 

With these properties in mind, we can now examine figure 4, which shows solutions for the 
three-fermion-line cases. (Again, all other three-fermion-line cases are related to the ones shown 
by parity or cyclic or reflection symmetries.) Using the properties of the vanishing vertices in 
figure |2j it is straightforward to show that only the desired QCD tree amplitudes on the left are 
produced by the gluino assignments on the right. The most subtle case is (3e). The two pairs 
of adjacent identical-helicity quarks force all the gluino flavors to be the same. However, this 
choice does not result in three distinct fermion lines in the pattern desired. There is a "wrong" 
fermion-line configuration which, fortunately, coincides with a permutation of case (3c). Hence 
we can subtract off this solution, as the second term on the right-hand side of (3e). 

We have also examined the cases with four different fermion lines. There are 16 cases, shown 
in figure 5, and each has a solution. Four of the cases, (4g), (4k), (41) and (4n), require two 
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Figure 5: The possible fermion-line configurations for amplitudes with four fermion 
lines. 

different terms, analogous to case (3e). Case (4p) requires three terms to remove all identical 
fermions; and case (4o) requires four separate terms. Case (4o) can be written as 

A QCB (^,qf,%,qt,qh% = A ^ ^"(l" 1+ > 1_ > 1+ > 1+ - ^ 1+ ' 1_ ) 

- [A x = 4 «™(2- 1+, r,2+, 2+, 2" 2+ 2-) - A M=i s™(2" 1+, 1", 2+, 2 + , 1", 1+, 2")] 

- [A*= A ^(2", 2+ 2-2+ 2+, 1", 1 + , 2-) - A*= 4 s™(2" 1+, 1~ , 2+, 2 + , 1", 1+, 2")] 

- A^= 4 SYM (2-, 1 + , 1", 2+, 2 + , 1", 1 + , 2-) (3.2) 
= A^= 4 syM (l", 1 + , 1", 1 + , 1+, I", 1 + , I") - A M = 4 SYM (2~, 1 + , 1-, 2+, 2 + , 2~, 2 + , 2-) 

- A^= 4 SYM (2-, 2 + , 2", 2+, 2 + , r, 1 + , 2~) + A M = 4 SYM ( 2 ~, 1 + , 1", 2+, 2+, 1", 1 + , 2~). (3.3) 

The first form of this equation indicates that three different wrong-fermion-line configurations 
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have to be removed. However, all removals can be accomplished with the help of different 
permutations of two other cases, (4f) and (41) respectively, 

A QCD (q^qt, ft, %,qt, ft, ft+ 9 2 + ) = ^ ^(1", 1 + , 2" 2" 1+ 1", 2+ 2 + ) , (3.4) 

^ QCD (ft, ?l" > 02, 97 > it, 97 94 > 92") 

= A^= 4 5yM (l", 1 + , 2~, 2", 2 + , 2-, 2+ 2+) - A M=i SYM (1~, 1 + , 2", 2", 1 + , 1", 2 + , 2 + ). (3.5) 

Case (41) itself requires a wrong-fermion-line subtraction. 

We have not yet ascertained whether any QCD tree amplitudes with more than four fermion 
lines are impossible to extract from M = 4 super Yang-Mills theory. Fortunately, for a fixed 
number of external partons, as one increases the number of fermion lines the number of Feynman 
diagrams decreases. Also, amplitudes with many external quarks typically contribute much less 
to a multi-jet cross section than do amplitudes with more gluons and fewer quarks. 

Finally, we remark on the conversion of pure-QCD tree amplitudes, that is amplitudes for 
quarks and gluons, into amplitudes that contain in addition a single electroweak vector boson, 
namely a W, Z or virtual photon. It is sufficient to compute the amplitude including the decay 
of the boson to a fermion-anti-fermion pair. 

Consider first the case of a virtual photon which is emitted from a quark q and decays to 
a charged lepton (Drell-Yan) pair £ + £~. We can extract this amplitude from a color-ordered 
amplitude with four consecutive fermions, as follows: 

A^{. . . , q + , I- , £+, g", . . .) = A QCD (. . . , q+, 9 2 ", ft+, 97, • • •) • (3.6) 

The color-ordering prevents gluons from being emitted from the quark line q 2 , or from the virtual 
gluon connecting qi and q 2 . Hence this virtual gluon is functionally identical to a virtual photon. 
The only other modification comes when dressing A 7 * with couplings. There is a relative factor 
of 2(—Q q )e 2 /g 2 when doing so, where the "2" is related to the standard normalizations of the 
QED interaction with coupling e, versus the QCD interaction with coupling g, and Q q is the 
electric charge of the quark q. (The lepton has charge —1.) 

It is possible to extract the amplitude (3.6) a second way, using one quark flavor instead of 
two, 

A 1 * (. • • , g + , r , £ + , q~ } . . .) = -A^ CD (. . . , g+ qf, ft, ft , • • •) • (3-7) 

This alternative works because the color-ordered interaction for g* — > q + q" is antisymmetric 
under exchange of q and q, and because the exchange of a gluon between identical-flavor quarks 
in the wrong channel is prevented by helicity conservation. 

If the virtual photon decays to other charged massless fermions, i.e. to a quark-anti-quark 
pair q'q', the only difference is of course to use the appropriate quark charge, —Q q e 2 — > Q q Q q e 2 . 
Because helicity amplitudes are used, it is also trivial to convert the virtual-photon amplitudes 
to ones for (parity-violating) electroweak boson production. For the case of combined exchange 
of virtual photon and Z boson, with decay to a charged lepton pair, the electric charge factor 
has to be replaced by 

2e 2 (-Q q + vi R vl R T z (s ei )), (3.8) 
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where v l L R are the left- and right-handed couplings of the lepton to the Z boson, v\ R are the 
corresponding quantities for the quark, 

p '« = 3 -MlUr z Mz < 3 ' 9 ' 

is the ratio of Z to 7* propagators, and Mz and Tz are the mass and width of the Z. 

Whether vl or v R is to be used in eq. (3.8) depends on the helicity assignment, i.e. on 
whether the Z couples to a left- or right-handed outgoing fermion (as opposed to anti-fermion) ; 
see ref. |37| for further details. The case of a W boson is even simpler, because there is no 
coupling to right-handed fermions (and no interference with another boson). 

The relevant MHV and NMHV amplitudes for four external fermions and the rest gluons, 
and for six external fermions and the rest gluons, have been converted in the above manner 
into tree amplitudes for Vqqg . . . g and VqqQQg . . . g, where V stands for W, Z or 7*. These 
NMHV amplitudes have been incorporated into the BlackHat library |31| and used there in 
conjunction with a numerical implementation |38| of the BCFW (on-shell) recursion relations |17| 
in order to obtain amplitudes at the NNMHV level and beyond. Including the MHV and NMHV 
formulae speeds up the numerical recursive algorithm by a factor of about four, in the present 
implementation. This approach was used to compute the real-radiative corrections entering the 
recent evaluation of pp — > W + 4 jets at NLO |27|, in particular the tree amplitudes for Wqq'ggggg 
and Wqq'QQggg. These amplitudes have nine external legs, after decaying the W boson to a 
lepton pair, so there are MHV, NMHV and NNMHV configurations, but no further. All seven- 
point configurations are either MHV or NMHV, so at most two recursive steps were required to 
hit an explicit formula (for example, in a schematic notation A 9 — > A 8 x A 3 — > A 7 x A 3 x A 3 ). 

We remark that the tree-level color-ordered amplitudes entering subleading-color loop ampli- 
tudes can have a more general color ordering from that required for purely tree-level applications. 
For example, in the pure QCD amplitudes with a single qq pair, only the color-ordered amplitudes 
in which the two fermions are adjacent are needed in eq. (3.1). However, the subleading-color 
terms in the one-loop amplitudes for qqg . . . g include many cases in which the two fermions are 
not color-adjacent, and the tree- level qqg . . . g amplitudes that enter their unitarity cuts have 
the same property. These color-ordered amplitudes are all available in Af = 4 super Yang-Mills 
theory, of course. 

Similarly, for computing subleading-color one-loop terms for single-vector boson production 
processes like Vqqg . . . g, one needs tree amplitudes such as A 1 * (. . . , q + ,g, £~, £ + , q~ ,...), in which 
the gluon g is color-ordered with respect to the quark pair, but not the lepton pair. These 
amplitudes are not equal to any particular color-ordered amplitude in M = 4 SYM, but one can 
generate them by summing over appropriate color orderings. For example, we have 

AT (..., q + ,g, r, £ + , q-, ...) = A^ CD (. . . , qt,g, qj, q+, ft, . • •) 

+ A QCD (...,g+, q^g,qt,^,...) + A <iCD (...,qt,^,qt,g,q^...). (3.10) 

The sum over the three permutations properly cancels out the unwanted poles as g becomes 
collinear with either i~ or £ + . 
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4 All gluon tree amplitudes 



In this section we present the general expression for an n-gluon tree amplitude, which we derive 
in section 7 from the solution of ref. |21 for a general M = 4 SYM super- amplitude. 

Amplitudes for n-gluon scattering are classified by the number of negative-helicity gluons 
occurring in them. Tree-amplitudes with fewer than two negative-helicity gluons vanish. In 
our conventions the gluon at position n is always of negative helicity, which does not present 
a restriction due to cyclicity of the color-ordered amplitude. The Parke- Taylor formula |11| for 
a maximally-helicity-violating (MHV) gluon amplitude, with the two negative-helicity gluons 
sitting at positions Co G [l,n — 1] and n, then reads 

4.(1+ • • • , (co - c -, (c + 1)+, . . . , (n - 1)+, n") := Af^o, n) = ^^^l) ' (41) 

with the total conserved momentum p = Y^=iPi- 

Next-to-maximally-helicity- violating amplitudes of degree p (N P MHV) then consist of (p + 2) 
negative-helicity gluons embedded in (n — p — 2) positive-helicity states. We parametrize the 
positions of the negative-helicity gluons in the ordered set (co, ci, . . . , c p , n) with q e [1, n — 1]. 

The general N P MHV tree-amplitude then takes the form 



^ PMHV (co, Cl ,...,c p ,n) = 


^ (4) (P) 

(12)(23)...(nl) 




x ^ 1 ' Rn;a\b\ 


r>{L2}]{U 2 } h{L P };{U p } /, „ pat h/ 
' iX n.;{/ 2 };a 2 62 ' • • • ' n n;{I p };apb p ' ^ UeX l C 0' " " 


\ 4 
• , Cp) ) 


all paths of length p 







(4.2) 



Let us now explain in turn the ingredients of this result, i.e. the sum over paths, the i?-functions, 
and the path- matrix E% . 

The sum over all paths refers to the rooted tree depicted in figure 6, introduced in ref. |21|. 
A path of length consists of the initial node "1". A path of length p leads to a sequence of 
p + 1 nodes visited starting with node "1". To clarify this all possible paths up to length p = 3 
are listed in figure 6. In general there are (2p)\/(p\(p + 1)!) different paths of length p, which 
is equal to the number of nodes appearing at level p in the rooted tree, since each final node 
unambiguously determines a path through the rooted tree. 

The i?-functions are generalizations of eq. (2.8) and may be written rather compactly with 
the help of eq. (2.6) as 

6 _L (q(a-l)) (6(6-1)) 

n,iIhab ' xlb (n {1} ba\a) (n {/} ba\a - 1) (n {/} a6|6) (n {/} a6|6 - 1) ' 1 ' 1 

they derive from the dual superconformal i?-invariants introduced in ref. |21|. In the above and 
in eq. (4.2), {/} is a multi-index deriving from the node in the associated path with the last 
pair of indices stripped off, e.g. {1%} = {61, ai, 6 2 , a 2 } for the last node of the first path of length 
p = 3. 
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Length p 


Paths 





[1] 


1 


[!]•[«! A] 


2 


[1] • [ai,&i] • [6i, ai; a 2 , 6 2 ] 
[1] • [ai,6i] • [a 2 ,6 2 ] 


3 


[1] • [ai, 6i] ■ [6i, ai; a 2 , 6 2 ] ■ [b u a u b 2 , a 2 ; o 3 , 6 3 ] 
[1] ■ [a 1; 6i] ■ [fei, ai; a 2 , 6 2 ] • [6i, ai; a 3 , 6 3 ] 
[1] • [ai, h] ■ [h, ai; a 2 , 62] • [a 3 , 63] 
[1] • [ai, 61] • [a 2 , b 2 ] ■ [b 2 , a 2 ; a 3 , 63] 
[1] ■ [ai,6i] ■ [02,62] • [03,63] 



Figure 6: The rooted tree encoding the sum over paths occurring in eq. (4.2). The 
table shows all paths up to length 3. 

In eq. (4.2) we used a further piece of notation, namely i?-functions with upper indices, which 
start to appear at the NNMHV level, and which we now define. Generally the i?-functions appear 
in the amplitude with an ordered summation over the last pair of indices, 

E Rn;{I h a b - (4-4) 

L<a<b<U 

i?-functions with upper indices indicate a special behavior for the boundary terms in this sum 
when a = L or b = U . Specifically if one has 

L<a<b<U 
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and the boundary of a summation is reached, then the occurrence of the spinor \a — 1) or \b) in 
the i?-function without upper indices (4.3) is replaced by a novel spinor depending on the upper 
indices as follows 



(L - 1| — y (£ L \ := (n\x nh x hh . . . xi p _ llp for a = L, 

(U\ — y (£u\ := {n\x nUl x UlU2 . . . x u Uq for b = U. 



(4.6) 



Effectively this amounts to the following formula for the upper-indexed 77-functions of eq. (4.5) 



■ryh,...,l p ;ui,...,Uq 
n n;{I};ab 



5 r , Kl) (n{I}ba\a-l) 
J^n;{I};ab ( n {/}fe a |^) (a(a-l)> 

6 , , (fe(b-l)) (n{/}afr|6) 
-n-n;{7};ab ( n {/} a 6|^) (6(6-1)) 



for a = L, 
for 6 = 17, 



D rn Kl) (n{/}6a|a-l) (&(&-!)) (n{J}a6|6) r _ r i , _ rj 

J^n;{I};ab (n{I}ba\£ L ) (a(a-l)> (n{/}a6|^> (6(6-1)) 1Ui " mlU U U ' 



(4.7) 



n;{7};ab 



else. 



In particular there is a term in the double sum where both a = L and b = U are reached 
and both replacements are to be made. The rules for constructing the sets of upper indices 
h, . . . , l p ] «i, . . . , u q in eq. (4.5) from the rooted tree are given in ref. |21|. 

To write down the path-matrix S^ ath we furthermore need to define the quantities 



( L= n)o6 
.' n '{bi,ai,...,b r ,a r };ab 



(rid) , 

{ribald) X[a,b-i](ci) - x 2 ab (na) X[m-i]( c *) 



(4.8) 
(4.9) 



(nMi . . . b r a r ab\ci) X[a,b-i](ci) - %lb ("Mi . . . b r a r |q) X[a r ,a-i](ci) , (4.10) 
with the support functions 

f 1 if i e fa, 61, , 



Now to every node [{7}; ab] along a given path and to every negative-helicity leg q we associate 
the entry of the path-matrix {S n )%y. ab . The entries (E n ) A B form a (p + 1) x (p + 1) matrix. 
Explicitly one has 



-n ath ( C 0; • • • ) c p) : = 



(nci) 



(ncp) 



\ 



7 ) c 
J «^ai6i 



,' r: «)a 1 i6i 



£7 ^ c p 



ai&i 



7 ^0 f£7 
Jn ^{/2};a 2 62 ^n;{/ 2 };a 2 6 2 



V 



' 77 ^0 
^ n /{7p};a p 6 p 



77 )C1 

^ n > {Ip};a p b p 



^ n nh};a 2 b 2 



Zn ) '{I p };a p b p J 



(4.12) 



Although R and S^ ath look rather involved at first sight, they are determined entirely through 
the external spinors Aj and Aj. 
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To clarify the construction principle let us write down the first three amplitudes in the N P MHV 
series explicitly: 



Ar HV (c , Cl ,n) 
^ 2MHV (co, Cl ,c 2 ,n) 



<12><23>... 


(nl> 






(12)(23)... 


(nl) 






(12)(23)... 


(nl) 



• (nc y 



2<a 1 <b 1 <n-l 



^ ^ Rn\a\b\' 
2<a\<b\<n 



(nc ) (nci) 



(4.13) 
(4.14) 



r: 



,ai+l<a 2 <fe2<fei 



(nc ) 



bl<a2,fe2<" 



'a 2 b 2 



> 


(nci) 




CO 
CO 

'bl,ai;a 2 6 2 


f- If 1 

V-^n/bi,ai;a2b2 


(nci) 


(nc 2 ) 


4 


V-^n/ a2 b2 


(^w)aifei 





(nc 2 ) 



(4.15) 



In appendix A we spell out the NMHV and N 2 MHV amplitudes explicitly. We provide a 
Mathematica package GGT with the arXiv.org submission of this article, which expands the 
master formula (4.2) explicitly for any choice of p, positions q and momenta AjAj. See appendix 
C for documentation. The formula (4.2) is implemented by the function GGTgluon. 

It should be mentioned that in practice the number of terms arising from the determinants of 
the path-matrix £;P ath is often quite small, see e.g. eq. (2.7). Moreover, for small n the number 
of non-zero terms in the nested sums can be relatively small. 



5 All single-flavor quark- ant i- quark- gluon trees 

Turning to the gauge-theory amplitudes involving massless single-flavor quark-anti-quark pairs 
we can write down a similarly general formula based on paths along the rooted tree of figure 6. 
In an abuse of notation, we refer here to a helicity +^ fermion as a quark, and a helicity — | 
fermion as an anti-quark. Looking at a color-ordered n-parton amplitude involving gluons and k 
quark-anti-quark pairs, g n ~ 2k (qq) k , it is again classified as a N P MHV amplitude by the number 
(2+p — k) of negative- helicity gluons. In such a color-ordered amplitude we furthermore consider 
an arbitrary ordering of the fermions. We then have 2+p+k 'special' legs (negative-helicity gluon, 
quark or anti-quark) in such an amplitude, whose position amongst the n legs we parametrize 
by the set 

(cq, . . . , c Ql , . . . , cgj . . . , c Qfc , . . . , cg fe , . . . , Cp + fc, n) . (5-1) 

Here q denotes the position of a negative-helicity gluon, c aj a quark and c§. an anti-quark 
location. Note that while the quark and anti-quark locations c ai and eg. are considered as ordered 
sets, i.e. c ai < c a . and eg. < cp. for i < j, there is no such ordering in the total set {c ai ,Cp.} 
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reflecting an arbitrary sequence of quarks and anti-quarks in the color-ordered amplitude. Again 
in our convention one negative-helicity gluon is always located on leg n 1 . 
The general N P MHV tree-amplitude for such a configuration then reads 



/± {qq) k ,n l C 0, • ■ • ) C «i ) • ■ • ) c ft ) 



" '"• , ' +A ""' 'l 2)(2 3)...(nl) 

3 



X 



X 



«W*0 E n^StS det(sr th |) det(«r h (^?) 



all paths \i=l 
of length p 



(5.2) 



Here sign(r) is the sign produced by permuting the quark and anti-quark legs into the alternating 
order {c ai , c^, c Q2 , c^j 2 , . . . , c afc , c^ fc }. 

Remarkably, the only difference from the pure gluon amplitudes is a change in the determinant 
factors of the path-matrix E% . With 2 + p + k 'special' legs the path-matrix associated to 
eq. (5.1) is now a (p + 1) X (p + 1 + k) matrix of the form 

/ (nco) 



;7path 



) 




• (raft) 


. (nc p+k ) 


CO 

<H&1 ' ' 


■ K— , n) ai i )1 


• ^"^0161 




CO 

{h};a 2 b2 


• ^ n ^{-r2}i026 2 • • 


• ^">'{/2};a26a • • 




CO 

{i"p};a p fe p 


• ^^{J p };a p &p • • 


■ y-n){I p };a p b p ■ ■ 


• ^J{J p };a p & p 



(5.3) 



The notation S^ 3 ^^ in eq. (5.2) now refers to the elimination of all the quark columns (the 
c aj 's) in the path-matrix and S^ ath (q g)|^ denotes the matrix with all the anti-quark columns 
removed (the Cg 's) after quark and anti-quark columns have been interchanged, i.e. cs, c ar 
The removal of k columns is of course necessary in order to turn the (p+l + k) x (p+1) matrix 
;rpath j n ^- S q Uare form, so that the determinant is defined. 

Let us again spell out some of the lower p amplitudes explicitly to clarify the formula (5.2): 



4 MHV/ \ 



^-qq,n ( C 0; c «i i c ^ i n ) 



(12)(23)...(nl) 
(12)(23)...(nl) 



(nc^) 3 • (nc ai ) 



(5.4) 



2<ai<fei<n-l 



/iNMHV/ 



(nc ) 
5^(p) 



(nco) (nc ai ) 



(5.5) 



(12)(23)...(nl) 



E 



2<ai<fei<n— 1 



*We comment in section 7 on how to circumvent this problem for the case without a single negative-helicity 
gluon. 
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'"Joi&i l^Jai&i 



(nc ai ) (nc a2 ) 



'oi&i 



A ^ HV (ca 1 ,c l ,c^ 1 ,c 2 ,n) 



(12)(23)...(nl) 



E 



n;ai&i 



2<ai<fei<n 



E 



^n\b\a\\a 2 b 2 



.ai+l<a 2 <b 2 <b 1 



(nci) 



^rc)ai&i 
'ST \ci 

v' 'rt/&i,ai;a2&2 



'0161 



(nc 2 ) 

(^n) 



'61,0150262 



C2 

Ol6l 

cr ^2 

^ n )b\,a\\a 2 b 2 



(ncx) 
(5 



(nc ai ) 



Cl 

ai&i 
— ' n ^ 01 ,a\;a 2 b 2 



Ca 1 
O161 



(nc 2 ) 



+ E 



_ 

0262 



6i<a2,fe2<« 



(nci) (nc h ) 



'61,0150262 
(nc 2 ) 



.^fOai&l 



■^ rft. 

"ft 



. 1= >0al&i 



a2&2 



^N 2 MHV^ 



" n )a2&2 
5< 4 )(p) 



C12&2 



"nJoi&l 
" ^2 

^~"nibi, 0150262 

(nci) 



JN MnV / _ _ \ ~ vr y 

(af) 2 ,n l c ai> c l> c /3i> c a 2 > c /3 2 > n J _ (12) (23) (nl) 



E 



R 



2<a\<b\<n 



E 



o0;ai6i 
0.5610150262 



_ai+l<a 2 <&2<6i 



(nci) 



^"Jal&i 
V'- J "/fc lj ai;a2&2 



+ E < 

6i<a2<62<'T' 



ai&i;0 
0262 



(nci) 

(s 

(■5 
(nci) 



«/ai6i 
n lb\, 0150262 



V" n /6l,ai;a262 

(nc ai ) 

V""/aibi 
x 1 ""^] 61,0150262 



(nc^) (nc^ 2 )_ 



^™)ai&i 
v ir: n)a 2 62 



\^n)a 2 b 2 



K'—n)a l bi 
\^n)a 2 b 2 



■ > — 1 \ Ca 

.""JoiSl 

w \ Cai 

^n)a 2 b 2 



n;aibi ' 



K^nia-Lbx 

\^—"fhJl ) y ,0150262 

(nc Q2 ) 



'-''Vai&i 
■~^ n /&i, 0150262 

(nci) 

(^"■)oi&i 
\^nJa 2 b 2 



K^nJarb! 
(W \ Ca l 
V" n /0262 



(5.6) 



(nc ai ) (nc 2 ) 



' 17 "\ c 2 

^™Jai&i 

' f7 "| c 2 
"W02&2 



(5.7) 



(nc ai ) (ric a2 ) 



«2 / 
/ £7 A Ca 2 



(5- 



This completes our examples. Some explicit formulae with the determinants expanded out may 
be found in appendix B, The formula (5.2) is implemented in GGT by the function GGTf ermionS. 
See appendix C for the documentation. 

Note that the master formula (4.2) reduces as it should to the pure-gluon scattering result (4.2) 
for a zero number of quark-anti-quark pairs, k — > 0. In that case no column removals are to be 
performed and the determinants combine to the power four. 
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6 All gluon-gluino tree amplitudes in Af = 4 super Yang- 
Mills 



The color-ordered tree amplitudes with fermions presented above were special in the sense that 
they apply both to massless QCD as well as N = 4 super Yang-Mills, due to the single-flavor 
choice which suppresses intermediate scalar exchange as argued in section 3. We now state 
the master formula for general gluino and gluon tree amplitudes in the M = 4 super Yang- 
Mills theory from which the above expressions arise. Through specific choices of external flavor 
configurations, however, this master formula may be used to produce color-ordered gluon and 
quark trees in massless QCD, as was discussed in section 3. 

Similar to the notation in section 5, for a general g n ~ 2k (qq) k amplitude with arbitrary flavor 
assignments to the 'quarks' 2 we have 2 + p + k 'special' legs (negative-helicity gluon, quark or 
anti-quark), whose position amongst the n legs we parametrize by the set 



(r r Al r? 1 

VKb • • • > ot\ i • • • i 



(6.1) 



Again the configuration of quarks and anti-quarks inside the gluon background may be arbitrary 
while the sets and are ordered. The general g n ~ 2k (qq) k amplitude with (2 + p — k) 
negative-helicity gluons is then expressed in terms of the i?-functions and the path-matrix H? 8 ^ 



defined above. It reads 



(<?<?) k ,n 



A™ V (co, 

(gg) K ,n \ u ' 



Bj v _ £ (4) (p) sign(r) 

,c h c p+k ,n)- {12) ( 23 y ^ nl) 



x 



all paths \i=l 
of length p 



4-fc 



£ sign(a) f[S A < B °u det ^P ath | (P*® <* <H 

(re5 fc i=l ^ Q 



(6.2) 



Here the notation S^ ath \ q refers to the path-matrix (5.3) with all 'quark' columns {c ai } removed, 
whereas 5P ath | ? (/% — > oij) denotes the same path-matrix where the 'anti-quark' column is 
replaced by one of the previously-removed 'quark' columns c av Also sign(r) is the sign of the 
permutation for bringing the initial color-ordering of the fermionic legs into the canonical order 
{c ai , cp l7 c a2 , Cp 2 , . . . , c ak , eg } and sign(cr) is the sign of the permutation of a. This formula is 
implemented in GGT by the function GGTf ermion. 



7 Proof of the master formula 

In this section we prove the master formula (6.2) for a general M = 4 super Yang-Mills tree 
amplitude with external gluons and gluinos of arbitrary flavor, as well as the more compact 

2 We refer to the gluinos in this and the following sections as 'quarks' in order to not introduce new terminology 
beyond that introduced already in sections 4 and 5. 
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single-flavor expression (5.2) and the pure gluon expression (4.2) as sub-cases. 

Amplitudes in M = 4 super Yang-Mills are very efficiently expressed in terms of a superwave 
function $ which collects all on-shell states of the PCT self-conjugate massless M = 4 multiplet, 
with the help of the Grassmann variables r\ A with A = 1,2,3,4 of the sit (4) R-symmetry, 

&(\, A, 77) =g+(X, A) + V A g A (X, A) + ^ V 4>ab{K A) + \^\ B rfe ABCD f{\ A) 

+ ^AV^f(A,A). (7.1) 

Here g ± are the ±1 helicity gluons, g A and the four flavor ±| helicity gluinos, and 0ab the 
six real helicity scalar states. The Grassmann variables r] carry helicity +^ so that the whole 
multiplet carries helicity +1. 

We can write the amplitudes in the on-shell superspace with coordinates (Aj, Xi,rji) |12,16,39|, 

An(\i, Aj, rji) = A($\... $ n ) . (7.2) 

Since the helicity of each supermultiplet <f>i is 1 the amplitude obeys 

hiAn{Xi, A f , 77*) = An{Xi, 7]i) , (7.3) 

where hi is the helicity operator for the ith leg. The component field amplitudes are then obtained 
by projecting upon the relevant terms in the rji expansion of the super-amplitude via 

9t ->V? = 0, 91 -> J d% = J dril drjf drf drf , ~g i)A -> J dr] A , gf ->■ - J d% rjf . 

(7.4) 

Note that the super-amplitude A n has a cyclic symmetry that can lead to many different but 
equivalent expressions in practice. In order to obtain compact expressions for component ampli- 
tudes often a judicious choice of this cyclic freedom can be made |21|. 

The general solution for tree super-amplitudes of Drummond and one of the present au- 
thors |21| takes the compact form 

jN p MHV _ £ (4) (p) £ (8) (g) ^ R {L 2 };{U 2 } R {L P };{U P } ^ 

(12U23) (nl) ' ' ^n,{/ 2 },a 2 6 2 ' " " " ' • fl n ) {/ P } ) o J ,bp ' - ) 

^ ' ^ / " " " \ 'all paths of length p 

where g aA = Y^=i ^tw is the total conserved fermionic momentum, and the dual superconfor- 
mal /^-invariant is 

Rn;{I};ab = Rn;{I};ab^ ^n;{I};ab(i) Vi^j j (7-6) 

in the notation of the previous sections. We now wish to project this result in on-shell superspace 
onto the relevant components for a general g n ~ 2k (qq) k amplitude. For this purpose we set all 
of the rji associated to positive-helicity gluon states to zero. This leaves us with the p + 2 + k 
remaining Grassmann numbers r\ c . associated to the 'special' legs of helicities —1 and ±|. To 
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project onto a negative-gluon state at position i one simply has to integrate eq. (7.5) against 
J d A r)i. Similarly, to project to a quark or anti-quark state at position i of flavor Ai one integrates 
eq. (7.5) against J drj i i or — J cftrjirjf* . All integrations have to be in color order. 

In accord with our convention above the leg n is chosen to be a negative-helicity gluon state, 
or an anti-quark if there are no negative-helicity gluons. This is a convenient choice because 
the only dependence of the super-amplitude on i] n is through the total fermionic momentum 
conserving delta function, which can be written as 



'p+k 
,i=0 



'p+k 



,i=l 



(go Cj) 

(cq n) 



'p+fc 

<* (4) (5>c,)ifc 

,i=0 



(7.7) 



For each path in eq. (7.5) the four-dimensional Grassmann delta functions in eq. (7.7), together 
with the p delta functions arising from the /^-invariants (7.6), may be written as 



p+i 



'p+k 



i=0 



0=0 
5 (4) 



'P+k , . 



'p+fc 



'p+fc 



i=l 



(c n) 



5 (4) 



i=o 



(7.8) 



,i=0 



i=l 



with the (p + 2) x (p + + 2) path-matrix . If we have a negative-helicity gluon at position 
n, the r] n integration is trivial and we can drop the trial r] n column and the row determined by 
the first delta function in eq. (7.8), ending up with the (p + 1) x (p + k + 1) path-matrix given in 
eq. (5.3). For the sake of readability we will drop the labels on the path-matrix in what follows, 
just denoting it by S, and assume a negative-helicity gluon at position n. The projection to the 
general g n ~ 2k (qq) k amplitude of eq. (4.2), with quarks of flavor A4 at positions c ai and anti-quarks 
of flavor Bj at positions Cp , 



1 U Cti J 



/i (gq) k ,n l C 0> ■ ■ ■ 

is then performed via the Grassmann integrals 



c p+ki n ) 



(-l) fc sign(r) 



p+k 

n 



J 



p+i 



'p+k 



k Vc- 
H /c /3( 



ip (4 me= 

j=0 



(7.9) 



i=0 



Here sign(r) compensates the minus signs we encountered by permuting the quark and anti-quark 
Grassmann variables from color order to the canonical order Yli=i I 



Let us first consider the pure gluon case 
r] c . — > E~^r\ c . immediately gives 

p _ 



i.e. k = 0. Performing the change of variables 




(7.10) 



i=0 
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thereby proving eq. (4.2). To evaluate the general integral (7.9) we first perform the (p + 1) 
four-dimensional integrations with respect to the 77's of the anti-quarks and gluons by making a 
change of variables similar to the pure-gluon case, leading to 

. k „ p+k k 

1=1 J i=o i=i 

... ,a fe } 

Here E\ refers to the elimination of all quark columns in the path- matrix. We can simplify the 
sum over % by making use of some basic facts of linear algebra. Namely, given a square matrix 
M = (rriij) with minors M^, its determinant and inverse can be written as 

det M = 1 )' • ',//,, det M i3 and (M~\ = ("l^'^f • (7-12) 

i 

Hence, eq. (7.11) simplifies to 

sign(r) (det ~ J] / <*< E det ^ a i)) < > t 7 " 13 ) 

i=i ^ i=i 

where ^1 (A — > ctj) denotes the replacement of an anti-quark column by a quark column. The 
remaining integrations are straightforward and give 

sign(r)(det^| g ) 4 ^ E sign(<r) J[ 8 AiB ^ det (s\ q {fi a{i) -+ a*)) • (7-14) 

o-e5 fc i=l 

The general M = 4 super Yang-Mills g n ~ 2k (qq) k amplitude is therefore 



(99) > n ^ t0 ' • • • ' c <* ' • • • ' % ■ ■ ■ ' W' n > (x 2) (2 3) . . . (n 1) 



T) 

X 



4-fe '•' 



e n ^su, ( det ^l) " e n ^ Bctw det ( s i,(a*> «o) • (7.i5) 



all paths \j=l / cGS^ i=l 

of length p 



Note that during the derivation of this formula we assumed that there is at least one negative- 
helicity gluon. The only change in the case k = p + 2 is that the r) n integration is no longer trivial 
and we put a gluino at position n. Hence, the path matrix has the size (p + 2) x [2p + 4) and 
eq. (7.15) still holds as its derivation did not depend on the matrix dimensions. Using (7.8) the 
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path matrix (4.12) then generalizes to the form 

(epci) (cqc 2 ) 







(c n) 



(con) 



{c c p+k } 
(c n) 



\ 





(nc ) 


(rid) 


(nc 2 ) 


. (nc p+k ) 







{- Y° 




(~ Y 2 


( \ c p+k 







(" Y° 

^ n '{I 2 };a 2 b 2 


{h};a 2 b 2 


(" Y 2 

\^ n '{I 2 };a 2 b 2 ■ 


( C Y p+k 







(~ Y° 




(~ Y 2 


( C7 \ c p+k 





(7.16) 



Generally the amplitudes take a more compact form if the gluino at position n is taken to be of 
helicity —1/2. Several explicit formulas for the MHV and NMHV cases can be found in Appendix 
B. In particular Appendix B.2.3 discusses a case without a negative helicity gluon at position n. 

As we are interested in QCD tree amplitudes, we need to decouple possible intermediate scalar 
states arising from the Yukawa couplings gA9B ( t )AB i n the AT = 4 super Yang-Mills Lagrangian. 
As discussed in section 3, one case in which this can be achieved (although not the only one 
needed for QCD) is when the external fermion states all have the same flavor, due to the anti- 
symmetry of 4> AB = e ABCD 4>cd- For this case, we specialize to Aj = Bi = A for all external 
fermion legs i in our master formula (7.15), and perform the sum over permutations explicitly, 
yielding 



.NPMHV/ \ 
/1 (gg) fe ,n l C 0> • • • j c «i j • • • ) c /3 j ■ ■ ■ > c p+k, n) 



8^>{p) sign( 



T 



(12)(2 3)...(n 1) 

N ' v ' v 'all paths \i=l 
of length p 

which reproduces eq. (5.2). 



D e n^a ; aA dets (^^L-( dets L) . ( 7 - 17 ) 
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A Explicit formulae for gluon trees 

Here we explicitly apply our formula (4.2) to the NMHV and NNMHV cases. 
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A.l NMHV amplitudes 

Without loss of generality, we take the negative-helicity gluons to be at positions Cq, Cx,n with 
c < C\. In the NMHV case only one path in figure 6 contributes and the path-matrix is a 2 x 2 
matrix whose determinant we denote by . Hence, the NMHV gluon amplitude is given by 



.4 



NMHV, 



c ) c l i n 



(12)...(nl) ^ 

N ' v ' 2<ai<6i<n-l 



R 



(d C0C1 V 



(A.l) 



where the determinant of the path-matrix is given explicitly by 



D 



ab 
n.st 



(n a) 



(n b) 



a<b 



' (n a)(nts\b) 
(n a)(b n)x 2 si 
(b a) (nts\n) 

K (n b)(nst\a) 



a < s <b < t . 
a < s <t < b . 
s < a,b < t , 
s < a <t < b. 



(A.2) 



For a > b one can use the antisymmetry of the determinant, D® b st = —D b ^ st . Equation (A.2) 
is exactly the result we already stated in eq. (2.7). This formula is implemented in GGT by 
GGTnmhvgluon. 



A.2 N 2 MHV amplitudes 

The negative-helicity gluons are taken to be a~ , b~, c~, n~ with a < b < c, without loss of 
generality. According to figure 6 there are two contributing paths. Denoting the determinants of 
their corresponding path-matrices by Df* and D^ bc , the NNMHV gluon amplitude is given by 



A 



N 2 MHV/ - l — „- - 



(a ,6 , c , n ) 



(12)...{nl)_.^ n;aifel ' 



2<ai<fei<n 



ai+l<a2<fe2<6l 



o0;ai6i 
n;b\a\\a-2,b-2. 



+ E <-X(A a6c ) 4 



bi<a,2<b2<n 



The explicit forms of the determinants of the path-matrices 

(n a) 



D 1 abc (n,a 1 ,b u a 2 ,b 2 ) :-- 



(n b) (n c) 

"n)axb\ ("n)a 1 6 1 ("rOaibj 

~ 'i/6i,ai;a2i>2 \"7i/&i,ai;a2&2 



J7 

.'—"n J [ bi ,ai;a2&2 



(A.3) 



(A.4) 



and 



Df c (n,a 1 ,b 1 ,a 2 ,b 2 ) :-- 



(n a) (n b) (n c) 

("n)ai&i ("Oaj&i (^fOai&j 
(^™)a 2 b2 ("™)<X2&2 ( 1 - 2 ")a 2 &2 



(A.5) 
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are given by 



abc 



{a n) (nbiai\b) (nb 1 aib2a 2 \c) 


a < ai < 


b,c <bi 


b < d 2 < c < b 2 


(n a)(nb 1 a 1 \b)(nb 1 a 1 \c)x 2 bo 


a < ai < 


b,c <bi 


b < a 2 , b 2 < c 


(a n)(6 c)(n6iaia2^2 ^&icn) 


a < ai < 


b,c <bi 


0-1 < 0-2 < b, c < 6 2 


(a n)(nb\ai \c) (nbiaia 2 b 2 \b) 

\ 1 \ - LJ -l/\ 1 / 


a < ai < 


b,c<bi 


di < d2 < b < b 2 < 


(n a)(c n)(nbiai \b)x 2 h x 2 nhn 


a < ai < 


b < 61 < 


c b < a 2 , b 2 < c 


(n a)(n c)x 2 r , h (nbid-id 2 bo\b) 


a < ai < 


6 < 61 < 


c d2 < b < b 2 


(a 6) (nfeiailn) (nb\d\b 2 d 2 c) 

\ / \ - LJ -| / \ _1_ J_ ^ | / 


d\ < 


a, b, 


c < 61 


b < d 2 < c < b 2 


(c b)(nb\ai \n) (nb\dib 2 d 2 la) 


a\ < 


a, b, 


c < 61 


d < a 2 < b, c < 62 


(6 a)(nbiai|n)(n&iai|c):r? 


a\ < 


a, b, 


c < 61 


b < d 2 ,b 2 < c 


(n6iai n) (x^ i, (nbia-i \a)(b c) 


d\ < 


a, b, 


c < 61 


d < d 2 < b < b 2 < c 


+ (nb^a^a2bo\b)(a c)) 










(a b)(nbidi \n)(nbia-\aobo\c) 


d\ < 


a, b, 


c < bi 


d 2 < d,b < b 2 < c 


(c 6)(n6iai n)(n6iaia2&2 a) 


ai < 


a, 6, 


c < 61 


di < d2 < a < &2 < & 


(6 c){nb\ai \a){nb\ai \n)x a2 b 2 


d\ < 


a, 6, 


c < 61 


0-1 < «2 < « < ^2 < & 


(c n) (a ftXnailnftOx^^x^^ 


di < 


a, 6 


< 61 < c 


6 < a 2 , &2 < &i 


(c n) (^ 262 (n6iai|a)(nai6i|6) 


di < 


a, b 


< h < c 


d < d 2 < b < b 2 


+ {na\b\ \a) {nbiaia 2 b2 \b)) 










(n c){a b){nbiaia2b 2 \na\bi) 


d\ < 


a, 6 


<h<c 


d 2 < d,b < b 2 


(c n) (nha^a) (nai&i \b)x 2 a2b2 


d\ < 


a, 6 


<h<c 


d < d 2 ,b 2 < b 


(n c) (naibi\b) {nbiaia 2 b 2 \a) 


d\ < 


a, 6 


<h<c 


d 2 < d < b 2 < b 



for 1 < di < a 2 < b 2 < bi < n and 

(n d) (nbidi\b) (nb2d2\c) 
(n d)(c n) (nhd^fyxl^ 
(n d)(b n) {nb 2 d 2 \c)x 2 aibi 
(n d)(b n)(n c)x 2 a ^x\^ 
{n d)(b c)(nb 2 d 2 \n)x 2 aibi 
,,,„. ) {na){cn){nb 2 a 2 \b)x 2 aibl 
(a 6)(nai6i|n)(n6 2 a2|c) 
(c n){a 6)(na 1 6 1 |n)x2 2b2 
(n b) (ndibi\d) (nb2d2\c) 
(n b)(cn)(ndibi\d)x 2 a2b2 
(c b) (nd\b\\d) (nb2d2\n) 
(n c) (nai&i | a) (na 2 &2 1 b) 



(A.6) 



D. 



d < d\ < b < b\ < c d2 < c < 6 2 
a<ai<6<6i<c &2<c 
a < ai, &i < 6 b < d 2 < c < b 2 
a < ai,bi < b b < d 2 ,b 2 < c 
a < ai,b± < b d2 < b,c < b 2 
a < ai,b± < b d2 < b < b 2 < c 
a 1 < d,b < bi < c a 2 < c < b 2 
a 1 < d,b < bi < c b 2 < c 
di < d < bi < b b < d 2 < c < b 2 
di < d < bi < b b < d2,b 2 < c 
di < d < bi < b a 2 < b, c < b 2 
di < d < bi < b a 2 < b < b 2 < c 



(A.7) 
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for 1 < ai < bi < 0,2 < 62 < n. For other orderings of a, b, c one can use the total antisymmetry of 
D^ bc and _Df fcc under permutations of a, b, c. It is quite astonishing that in 28 out of 30 cases these 
determinants are given by a single term. This formula is implemented in GGT by GGTnnmhvgluon. 

B Explicit formulae for trees with fermions 

Here we explicitly write out our formulas (5.2) and (6.2) for the MHV, NMHV and NNMHV 
cases with up to six fermions. 

B.l MHV amplitudes 

The simplest amplitudes involving fermions are the MHV amplitudes. The amplitudes with one 
negative-helicity gluon a and two fermions of opposite helicity and the same flavor, b (+|) and 
c (— |) (or vice- versa), are given by 

Ma-Ac)=^p) { ^ ?{a ^ iy (B.l) 
Ma-M^-S^)^^. (B.2) 

These formulae correspond to case (1) in figure 3. We note that the latter formula is related to 
the former one by a reflection symmetry, under which the cyclic ordering is reversed and there 
is a relabeling b c. In the NMHV case we will omit formulae that can be obtained from the 
presented formulae by a reflection symmetry. 

An equally compact formula can be obtained for the MHV amplitudes with four fermions and 
only positive-helicity gluons: 

A n (a A ,b B ,c c ,d D ) = * ( y ){b , d) * {5 AB 6 CD (da}(cb)-6 AD 6 BC (dc}(ab}) , (B.3) 

(1 2) ... (n 1) 

AJa A ,b B ,c c ,d D ) = .fy^, (5 AD S BC (d b)(a c) - 5 AC S BD (d a)(b c» , (B.4) 

(1 2) ... {n 1) 



which in the single-flavor case simplifies to 



A n (a, b, c, d) = (1 ^ X w x , (B.5) 

- 5^(p)(cd) 3 (a b) 
A n (a,b,c,d) = j^-j . (B.6) 



Equation (B.6) corresponds to case (2a) in figure 3, whereas eq. (B.3) for A = B ^ C = D 
corresponds to case (2b). 
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To complete the list of MHV amplitudes with up to four fermions we also give the MHV 
amplitude with four positive-helicity fermions and one negative-helicity gluon: 

A n (a A ,b B ,c c ,d D ,n~) = [ d V A [ I dtf J d^J d^ n AT V 

n a) (n b) (n c) (n d) . (B.7) 



S^{p)e ABCD 



(12). ..(n 1) 
This amplitude is not needed for QCD. 

B.2 NMHV amplitudes 
B.2.1 Two fermions 

To illustrate the use of our master formula (5.2) we compute the NMHV amplitude with two 
opposite-helicity fermions at positions a, b and two negative-helicity gluons at positions c and n. 
At this stage we leave the color order arbitrary. Starting with the path-matrix 



;7path 



(B.8) 



(n c) (n a) (n b) 

("n)st {^njtt ("n)rf 

we can immediately write down the amplitude 

(a nNMHV_ £ (4) (p)sign(r) , y 

{ A n) q q ~ Q 2 ) . . . (n 1) ^ ^^st [^st ) , 



<s<t<n 



where the 2x2 determinant D® b st has been defined in eq. (A.2). As already stated, the last 
equation holds for an arbitrary color ordering. In the following we take a < b < c and specify 
the color ordering: 



At -A J(4) (P) 

A n (a,b ,c,n 



(1 2)...(nl)L 



ab)(bc) 3 ^2 (nts\n) Rn,st 

Ks<a,b,c<t<n 

-(bc) 3 (an) (nts\b}(nts\n} 3 R n , st 

a<s<b,c<t<n 



— (c n) 3 (a n) 

Y (nst\bf(nts\b)R n 

a<s<b<t<c 

— (c n) 3 (a b) Y (nst\b) 3 (nts\n)R nt 



,st 

l<s<a,b<t<c 



(B.10) 



A n (a,b,c ,n ) 



(12}...(n 1)L 



+ (ac)(bc) 3 Yl (nts\n) A R n . 



St 



Ks<a,b,c<t<n 
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+ (an)(bn} 3 ^ (nts\c) 4 R nM 

b<s<c<t<n 

+ (nc) 4 (an)(bn) 3 ]T (x 2 t ) 4 R n>st 

b<s<t<c 

+ (cn) 4 (nst\b) 3 (nst\a)R n)St 

l<s<a,b<t<c 

+ (bc} 3 (an) Y {nts\c)(nts\n) 3 R n 



St 

a<s<b,c<t<n 



+ (cn) 4 (an) ^ x 2 st (nst\b) 3 R nyS 



a<s<b<t<c 



(B.11) 



These simplified expressions are implemented in GGT by GGTnmhv2f erm (see appendix C). 
B.2.2 Four Fermions 

We proceed with the NMHV amplitude with four fermions at positions af 1 , a Ai , bf 1 , b B 2 and one 
negative-helicity gluon. Without loss of generality we put the negative-helicity gluon at position 
n. Again we leave the color ordering arbitrary. A straightforward application of our formulas 
(6.2) and (5.2) yields 



(a \ NMHV _ ^ •'(P)sign(r) ^ ( r) blt > 2 \ 2 ( R A \ B \R A i B 2 n aib ' 2 n bia ^ £A 1 b 2 jtA 2 b 1 n a 2 b 2 n b iai 

v n Hw) 2 ~ (12) (n I) " ' ,St V n ' st J \ n > st n ' st u n\st u n;st 

(B.12) 

in the Af = 4 super Yang-Mills case, and 

(a nNMHV_ ^ 4 )(p)sign(r) ^ ~ f n a ia2 /r> i o\ 

\ / ' ' ' \ / ± <s <t<n 



for single-flavor QCD, with D®. st defined in equation (A.2). Taking a < b < c < dwe now specify 
the color ordering, 

A n (a A ,b B ,c C 4 D ,n-) = f )( P } x 

(1 2}...(nl) 



x 



+ 5 AC 5 BD (an)(bc)(cd) 2 ^ (nts|d) (ntsln) 3 ^ - (c ^ d) 

a<s<6,d<t<n 



+ 5 AC (5 BD (ac/)(6c)(rfc) 2 (^(n) 4 ^ - (a b) 

Ks<a,d<t<n 

+ 5 AC 5 BD (na}(bc)(dn} 3 x 2 st (nst\n) (nst\c) 2 R n , st 



a<s<b,c<t<c 
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-5 AD 5 BC (an)(dn) 3 ^ (nts\c) (nst\b) (nst\c) 2 R HtSt 



a<s<b,c<t<d 

+ 5 AC 5 BD (bc)(dn) 3 (nst\a)(nts\n)(nst\c) 2 R ntSt - (a^b) 

Ks<a,c<t<d 

(B.14) 

where "(c -H- d)" implies the substitution c -H- d in the arguments of the spinor strings, as well 
as the corresponding substitution C -H- Z? in the arguments of the 5 functions, but no change in 
the summation range. The other inequivalent orderings of quarks and anti-quarks are, 

A n (a A ,b B ,c c ,d D ,n-) = f )( P } x 

(12)...(ral) 



x 



+ S AB S CD {an){bn) 3 ^ (nts\c){nts\d) 3 R n , st 

b<s<c,d<t<n 

+ 5 AB 5 CD {n a) (c b) (b d) 2 ( nts \ d ) (nts\n) 3 R HySt - (b <* d) 



a<s<b,d<t<n 



+ 5 AB 5 CD (ad)(bc)(db} 2 (nts\n) A R niSt -(a^c) 

Ks<a,d<t<n 



+ 5 AB 5 CD (na){dn) 3 (nb) 3 {nts\c){x 2 st ) 3 R n , st 

b<s<c<t<d 

+ 5 AB 5 CD (na}(nc}(nb} 3 (nd} 3 ]T (x 2 st ) 4 R n , st 

b<s<t<c 

+ 5 AB 5 CD (na}(bc)(dn} 3 x 2 st (nst\n}{nst\b} 2 R^ st 

a<s<b,c<t<d 

-S AD 5 CB {na){dn) 3 ^ (nts\b) (nst\c) (nst\b) 2 R nM 

a<s<b,c<t<d 

+ 5 AB 5 CD {na){nc){dn) 3 ^ x 2 st (nst\b) 3 R r< 



L n,st 

a<s<b<t<c 



+ 5 AB 5 CD (bc)(dn) 3 (nst\a)(nts\n)(nst\b) 2 R ri:St - (a <-> c) 

Ks<a,c<t<d 



+ 5 AB 5 CD (nc)(nd) 3 ^ (nst\a) (nst\b) 3 R n>st 



l<s<a,b<t<c 



(B.15) 



A n (a A ,b B ,<?,d D ,n-) = - .* A)(P ) x 

(1 2) ... (n 1) 
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X 



+ 5 AB 5 CD (an)(bn) 3 (nts\d)(nts\c) 3 R ntSt 

b<s<c,d<t<n 

+ 5 AB 5 CD (n a){n d}(b n) 3 x 2 st {nts\c) 3 R nM 

b<s<c<t<d 

+ 5 AB 5 CD (na){db)(bc) 2 (nts\c)(nts\n) 3 R nM - (b <->■ c) 

a<s<b,d<t<n 

+ 5 AB 5 CD (na)(nd)(bc) 2 ^ (nts\c) (nst\b) (nts\n) 2 R n , st - (b <-> c) 

a<s<b,c<t<d 

+ 5 AB 5 CD (ac}(bd}{cb} 2 (nts\n} 4 R nySt -(a^d) 



Ks<a,d<t<n 



+ 5 AB 5 CD (ac)(nd)(cb) 2 (nst\b) (nst\n) 3 R nM - (b ^ c) 

Ks<a,c<t<d 



+ 5 AB 5 CD (na)(nd)(nb) 3 (nc) 3 £ (x 2 st ) 4 R n , st 

b<s<t<c 

+ 5 AB 5 CD (na}(nd}(cn} 3 x lt( nst \b) 3 Rn,st 

a<s<b<t<c 

+ 5 AB 5 CD (nd)(nc) 3 (nst\a) (nst\b) 3 R HySt 



Ks<a,b<t<c 



(B.16) 



A n (a A ,b B ,c c ,d D ,n-) = - h ^ ){p) , n x 

(12)...(nl) 



x 



+ 5 AB 5 CD (bn)(an) 3 ^ (nts\c)(nts\d) 3 R^ st 

b<s<c,d<t<n 

+ 5 AB 5 CD (bd)(an) 3 (nts\n)(nts\c)(nts\d) 2 R^ st - (b c) 

a<s<fe,d<t<n 

+ 5n CD (y}(ac}(cfa) 2 ("*s|«) 4 £n,«f - (& «-> c) 

Ks<a,a!<t<n 

+ <W™(n 6><d n) 3 (n a) 3 £ (n^c)^) 3 ^ , si 

b<s<c<t<d 

+ 5 AB 5 CD (nb)(nc)(na) 3 (nd) 3 £ (a^) 4 ^ 



b<s<t<c 



+ 5 AB 5 CD (na) 3 (nd) 3 ^ (^^(n^lc)^) 2 ^ - (6 O c) 

a<s<6,e<t<<i 
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+ 5 AB 5 CD (nc)(dn) 3 (na) 3 ^ {nst\b){x 2 st f R ntSt 



a<s<b<t<c 



+ 5 AB 5 CD (ac)(dn) 3 ^ (nst\b) (nts\n) (nst\a) 2 R n>st - (b <* c) 

Ks<a,c<t<d 

+ 5 AB 5 CD (nc)(nd) 3 ^ (nst\b) (nst\a) 3 R nM 

l<s<a,b<t<c 



(B.17) 



For A ^ B, and all fermions cyclically adjacent we have 
rB , _ D , —a 5( 4 )(p) 



A„(a A ,(a + l) , (a+2)", (a+3) ,n") = -- 



x 



(12)...(nl) 

(a n)(a+2 a+3)(a + l a+3) 2 (nta+2|a + l)(nta+l|n) 3 J R niftm 



a-h3<t<n 



+ (a+2 a+3)(a a+l)(a+l a+3} 2 (nts|n) 4 i? niSi 

Ks<o,st3<«n (B.18) 

+ (a n)(a+3 n) 3 (a+l a+2) 4 (n|x nttf3 |a+2](n|a;„^„i|a + l](n|a; na+ i|a+2] 2 ^ nia+ i a+ 3 

— (a a+1) (a+3 n) 3 (nsa+3|a+2)(nsa+3|n)(nsa+3|a + l) 2 _R niSQ+3 
ks<o 

This amplitude may be used to generate the NMHV amplitudes for Vqqg . . .g, as discussed in 
section 3. 

In the single-flavor case we obtain 



- 5 (4) (p) 
A n (a,b,c,d,n ) 



(1 2) ... (n 1) 



+ (na)(cd) 3 ^ (nts \ b) (nts \ n) 3 R n , s t 

a<s<b,d<t<n 

+ (ab)(dc) 3 J2 (nts\n) A R ntSt 

Ks<a,d<t<n 

+ (na)(dn) 3 ^ (nts\b) (nst\c) 3 R n;St 



a<s<b,c<t<d 



+ (ba)(dn) 3 ^ (nts\n) (nst\c) 3 R n>st 



Ks<a,c<t<d 



(B.19) 



- _ 5 {4) (p) 
A n (a,b,c,d,n ) - 



(1 2)...<nl) 



+ (an)(bn) 3 ^ (nts\c)(nts\d) 3 R nM 

b<s<c,d<t<n 

+ (na)(db) 3 ^ (nts\c) (nts\n) 3 R n>st 



a<s<b,d<t<n 
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+ (ac)(bd) 3 (nts\n) 4 R nM 

Ks<a,d<t<n 

+ (na}(dn} 3 (nb} 3 (nts\c}(x 2 st ) 3 R n , st 

b<s<c<t<d 

+ (na)(nc)(nb) 3 (nd) 3 ]T (x 2 st ) 4 R n , st 

b<s<t<c 

+ (na)(nd} 3 ^ {nts\c)(nst\b) 3 R nM 

a<s<b,c<t<d 

+ (na)(nc)(dn) 3 x 2 st (nst\b) 3 R nM 

a<s<b<t<c 

+ (ac)(dn) 3 ^2 {nts\n) (nst\b) 3 R nM 

l<s<a,c<t<d 

+ (nc)(nd) 3 ^ (nst\a)(nst\b) 3 R n , st 



Ks<a,b<t<c 



<J (4) (p) 

A n (a,b,c,d,n ) = 



(1 2) ... (n 1) 



+ (an)(bn) 3 ^ (nts|d) (nts|c) 3 i? njSt 

&<s<c,c(<t<n 



+ (n a)(n d)(b n) 3 ^ x 2 st (nts\c} 3 R n:St 

b<s<c<t<d 

+ (na)(cb} 3 ^ (nts\d) (nts\n) 3 R nM 

a<s<b,d<t<n 

+ (n a)(n d)(b c) 3 ^ x 2 t (nts\n) 3 R n>st 



a<s<b,c<t<d 



+ (ad)(bc) 3 (nts\n) 4 R nM 

Ks<a,d<t<n 

+ (b c) 3 (n d) ^ ( nst I a ) ( nst I n ) 3 Rn,st 

Ks<a,c<t<d 



+ (na)(nd)(nb) 3 (nc) 3 ^ (4) 4 J R„, st 

b<s<t<c 

+ (na}(nd}(cn) 3 ^ x^ns^fe) 3 ^^ 



a<s<b<t<c 



(B.20) 
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A n (a,b,c,d,n ) 



(1 2) ... (n 1) 



+ (nd)(nc} 3 ^ (nst\a}(nst\b} 3 R nt 

l<s<a,b<t<c 

+ (bn)(an) 3 V) (nts|c) (nts\d) 3 R n 

b<s<c,d<t<n 



st 



st 
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+ (bc)(an) 3 (nts\n)(nts\d) 3 R n)S t 

a<s<b,d<t<n 

+ (bc)(ad) 3 (nts\n) 4 R 

n,st 

Ks<a,d<t<n 

+ (n b)(d n) 3 (n a) 3 ^ (nts\c){x 2 st ) 3 R n , st 

b<s<c<t<d 

+ (nb)(nc)(na) 3 (nd) 3 ^ (o4) 4 #n,st 

b<s<t<c 

+ (na} 3 (nd} 3 (bc) (nts\n)(x 2 st ) 3 R n , st 

a<s<b,c<t<d 

+ (n c)(d n) 3 (n a} 3 ^ (nst\b)(x 2 st ) 3 'R n , s t 

a<s<b<t<c 

+ (bc)(dn) 3 S2 (nts\n)(nst\a) 3 R nySt 

Ks<a,c<t<d 

+ (nc)(nd) 3 {nst\b){nst\a) 3 R nM . (B.22) 

Ks<a,kKc 

These simplified expressions are implemented in GGT by GGTnmhv4f ermS for the single-flavor case, 
and by GGTnmhv4f erm for the general-flavor case. See appendix C for the documentation. 

B.2.3 Six fermions 

In the case of the six-fermion NMHV amplitude there is no negative-helicity gluon for us to put at 
position n as we did in the previous examples. The fermions are at positions af 1 , a^ 2 , a^ 3 , 5f 1 , &f 2 
and n B ' A . This time the path-matrix (7.16) is given by 



;7path 



■y (&_2 ai) (ft_2 02) (b_2 gg) ^ 
{62 n) (b2 n) (b'2 n) 



( & 5l > 

(62 n) 

(n b\) (n 6 2 ) (n a x ) (n a 2 ) (n a 3 ) 

\ ( 1= n)st (^rOst ("n) s t {^njst {^rijst J 

As ingredients of formula (6.2) we need the determinants 

det (~P at \) = , det {S^Xih a,)) = D% 



0,62 



(B.23) 



(B.24) 
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det [S^ th \ q (b 2 at)) = D 



bioj 
n:st 5 



det (S path \ 



in — > a 



0) 



We recall that -D^: st has been defined in eq. (A. 2) and the 3x3 determinant D°^ c st reads 



j~\abc 

U n\st 



(M 

(bn) 

(n a) 







(n b) 



(b c) 
(bn) 

(n c) 



For a < b < c we have 



a 6)(nts|c) 



(a &>(~ n )^ + (& c)(~ n )« +(c a)(~ n 



'si 
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U n\st 



St 



a b){c n)x 
nts\a) (c b) 
n a)(b c)x 
a b)(c n)x 2 st — (a c)(nts\b) 
a b) (nst\c) 
c b)(nst\a) 



b < s < c < t 
b < s < t < c 
a < s < b,c < t 
a < s <t <b 
a < s < b < t < c 
s < a,b < t < c 
s < a < t < b , 
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and D^ t is totally antisymmetric in a, b, c. Thus, the M = 4 super Yang-Mills NMHV six-fermion 
amplitude is 



/ a nNMHV 



sign(r) 



(12>...(nl> 1 



E a 

<s<t<n 



p\bib 2 

n;stJ-^ n:st 



^n;si -L^mst ^njst 

_ A A 1 S 3A A 2 B 2A A 3 S 1 n a 3&2 n^i a 2 n&ifeai 
iy n;s j -t^n;st ^ n:st 
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^A 1 S 1 ^A 2 B 3 ^A 3 B 2 £) a l b 2 £) b l a 3 £) b lb 2 a 2 



+ S A lB3S A 2 B lS A 3 B 2 D a 2 b 2 D b,aj pl&ai 



which in the single-flavor case (5.2) reduces to 

vNmhv 5 (4) (p)sign(r) 



(An) 



(qq) 3 



(12)...<nl) 



E * 



t-)6i6 2 



nii«2<i3 



(B.29) 



L<s<t<n 

These simplified expressions are implemented in GGT by the functions GGTnmhv6f ermS for the 
single-flavor case and GGTnmhv6f erm for the general-flavor case. See appendix C for the docu- 
mentation. 



B.3 N 2 MHV amplitudes 
B.3.1 Two fermions 

We continue the list of quark-gluon amplitudes by applying the master formulas (5.2) and (6.2) 
in the N 2 MHV case with up to six fermions. The amplitude with three negative-helicity gluons 
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at positions ci, C2, n, a quark at position a and an anti-quark at position b, is 
8^(p) sign(r 



/ a nN 2 MHV 



(1 2) ... (n 1) 



2<a\<b\<n 



. ai+l<a 2 <62<fei 

+ E *5£2 ■ A ciC2a (a ciC25 ^ 



cic 2 b 



6l<a2<62<« 

with the 3x3 determinants and D^ bc from eqs. (A. 6) and (A. 7). 
B.3.2 Four fermions 



(B.30) 



For the amplitude with two negative-helicity gluons at positions well as quarks and 

anti-quarks at positions af 1 , a^ 2 and /Jf 1 , /3-f 2 , we obtain 



, sN 2 mhv _ £ (4) Q) sign(r) ^ ~ 
1^(^)2 - /i2)...(nl) ^ ^n;oi6i 



X 



2<ai<bi<n 



E (A^J (^SA^A* 08 - ^A^A*" 1 ) (B.31) 

ai<a 2 <6 2 <6i 



I \ D a l fe i;0 / T-)C^i/3 2 \ / XB 1 rB 2 7-jcai/3 2 t-)C/3io 2 _ rs 2 rSj n co 2 /3 2 ncftai 

^ -"-7150262 1^2 I l ^! ^ 1 ^ Mi ^1 A 2 "j2 -^2 

6i<a 2 <6 2 <n 



in the M = 4 super Yang-Mills case, and 

(r) v 5 

j\ / , , - n *ra;ai&i 



(a xN 2 mhv _ <5 (4) (p) sign(r) 
1 n) ^ 2 " (12)...(nl) 



2<«i <6i <n 



E ^SS^ • A CQ1Q2 (A* A ) (B.32) 

ai<a 2 <6 2 <6i 



+ E k tl ■ A ClC2a (A C1C25 

6l<<!2<62<™ 



for single-flavor QCD. 



B.3.3 Six fermions 



For the M = 4 super Yang-Mills amplitude with one negative-helicity gluon at position n, quarks 
and anti-quarks at positions af 1 , a^ 2 , af 3 and /3f \ ft® 2 -, P3 3 , our master formula yields 



(A. 



>N 2 MHV 



S^(p) sign! 
(1 2>...(nl> 



-j\ / - n -7i5ai6i 



X 



2<ai<6i<n 



35 



X 



E *5&;-h ( *S*2*2 A^ 2 A^ 2 A^* 3 ± permutations of 

ai<a2<b2<bi 

+ J2 D? M [5 B A \ 5 B A \ 5*1 D« lhh D^ a2h D^ aa ± permutations of 



OLi 



b\<a,2<b2<n 



OLi 



(B.33) 

which in the single-flavor case simplifies to 
(12)...(nl) 



/ , xN 2 mhv _ £ (4) (p)sign(r) ^ _ 



rl n;6iai;a 2 62 / 

2<ai<6i<n L ai<a2<b2<?>i 

Epai6i;0 naia2a3 / 7~)/3l " X 

^50362 ^2 ^2 



-" , n;a 2 b2 ^2 I ^2 

bl<a2<&2<« 

(B.34) 

We recall that these formulas hold for arbitrary color-orderings of the n partons. 

We have implemented all of the above simplified expressions for NNMHV amplitudes with 
up to six fermions in the functions GGTnnmhv2f erm, GGTnnmhv4f erm, GGTnnmhv6f erm in the GGT 
package. 



C The Mathematica package GGT 

Here we describe the Mathematica package GGT (gluon-gluino trees) provided with the arXiv. org 
submission of the present paper and also accessible via http://qft.physik.hu-berlin.de. 

The idea is to provide the formulas derived in the present paper in computer-readable form, 
such that the interested reader can use them without having to type them in. We have also 
included a simple numerical evaluation routine for given phase-space points in the GGT package, 
as well as an interface to the spinor-helicity package S@M |42|. The issue of computer speed opti- 
mization will be commented upon below. 

Let us now describe the different functions in GGT and then give a specific example. The 
following functions are provided in GGT 

• GGTgluon [n,H] 

gives the n-gluon amplitude (4.2), with the positions of the negative-helicity gluons given 
by the list H. 

• GGTf ermionS [n, gluonlist, fermlist, afermlist] 

gives the n-parton amplitude (5.2) of an arbitrary number of gluons and single-flavor 
fermion/antifermions. The positions of the negative-helicity gluons, helicity +| fermions, 
and helicity — | anti-fermions are given by the lists gluonlist, fermlist, and afermlist, 
respectively. 



36 



• GGTf ermion [n, gluonlist, fermlist, afermlist] 

is the generalization of GGTf ermionS to multiple fermion flavors, eq. (6.2). The positions 
of the negative-helicity gluons are given by the list gluonlist. The positions qi,q~i and 
flavors Ai,Bi of the helicity +| fermions and helicity — | anti-fermions are given by the 
lists f ermlist= {{qi, A t }, . . .}, and af ermlist= {{q~i, Aj}, . . .}, respectively. 

• GGTsuperamp [n, A;] 

is the N fc MHV superamplitude of n superfields, with the MHV superamplitude factored 
out, in terms of the R invariants. 

Let us give an example. We can load the GGT package using 

« GGT.m 

Suppose we want to evaluate a gluon amplitude. Typing 
GGTgluon[6,{3,5,6}] 

prints the 6-gluon NMHV amplitude with helicity configuration -H 1 , 

1 

(1|2)(2|3)(3|4)(4|5)(5|6)(6|1) 

/ (2|1)(4|3) (s 2|4 (6|3)(6|5) + (6|x M |x 4 , 2 |3) (6|5)) 4 
\S2,4 (6|x 6i2 |x 2 ,4|3) <6|ar 6j2 |aJ2,4|4) (6|x 6i4 |x4 i2 |1) (6|x 6)4 |z 4i2 |2) 
| (2|1)(5|4) (s 2 , 5 (6|3)(6|5) (6|x 6|5 |x 5 , 2 |3) (6|5)) 4 

S2,5 (6 1 £6,2^2,5 1 4) <6|^ 6;2 |^2,5|5) (6|x 6 , 5 |x 5i2 1 1) (6 |x 6j5 1 a; 5; 2 1 2) 
| (3|2)(5|4) (s 3i5 (6|3)(6|5) + (6|x 6 , 5 |x 5 , 3 |3) (6|5)) 4 \ 

s 3) 5 (6|x 6>3 |x 3i5 |4) (6|x 6i3 |x 3: 5|5) {6\x 6t5 \x 5 , 3 \2) (6|x 6)5 |rr 5:3 |3) / 

GGT formatted the output for better readability. The underlying formula, which can be ac- 
cessed explicitly, e.g. by using Inputf orm[. . .] , depends on the following quantities: The 
spinor products (ij) are denoted by GGTspaa[i, j] . Differences between dual coordinates Xij = 
Pi + Pi+i + • • • + Pj-i are denoted by GGTx[i, j] . Finally, the abbreviation xf^ = Sjj-i is used 
and denoted by GGTs[i,j-l]. 

In order to obtain numerical values, we can use the spinor-helicity package S@M |42|. The 
function GGTtoSpinors converts the expression into one that can be evaluated by the latter 
package. In our example, the commands 

« Spinors.m 
GenMomenta [1 ,2,3,4,5,6] 

load the S@M package and use one of its functions to generate arbitrary momenta for a six-particle 
scattering process. Finally, numerical values of the amplitude at that phase-space point can be 
obtained by the command 

GGTtoSpinors [GGTgluon [6 , {3 , 5 , 6}] ] //N 
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A faster implementation for the numerical evaluation of the GGT formulas is provided by the 
function GGTgenvar [P] which generates the spinors and region momenta for a numerical evalu- 
ation of an amplitude at a desired phase-space point P = {p%,P2, ■ ■ ■ ,p n }- F° r example, for the 
kinematic point given in eq. (4.6) of ref. |43| (which to save space we give here to only three 
significant digits), one would use 

GGTgenvar [{{-3.0, 2.12, 1.06, 1.84} , {-3.0, -2.12, -1.06, -1.84} , {2.0, 2.0, 0.0, 0.0} , 
{0.857, -0.316, 0.797, 0.0} , {1.0, -0.184, 0.465, 0.866} , {2.14, -1.5, -1.26, -0.866}}] 

One can then evaluate an amplitude numerically by the command 

GGTnumeric [GGTf ermionS [6, {1, 6}, {2, 4}, {3, 5}] ] 
- 0.496838 + 0.0714737 i 

This approach is considerably faster than the GGTtoSpinors [. . .] //N function discussed above. 

Let us comment about the evaluation time needed using our approach. It is clear that for any 
serious applications or for comparisons with other methods, one should implement our analytical 
formulas using a low-level programming language, such as C, C++ or FORTRAN. For example, 
an implementation of the NMHV formulas in C++ results in a speedup of orders of magnitude 
over a similar implementation in Mathematica. Moreover, it is important to efficiently cache 
(store the numerical values of) quantities that are used repeatedly. In this spirit, the Mathematica 
demonstration package GGT provides a computer-readable version of the formulas needed for such 
an approach, so that the user does not have to type them in manually. 

Our analytical formulas are very similar, and in some cases identical, to the ones obtained 
in a very recent paper |44|. The latter also correspond to solutions of the BCFW recursion 
relations, based on refs. |45|, but may differ in form since they can correspond to different 
factorization channels. Another difference is that they are written using momentum-twistor 
variables |46|. Ref. |44| contains a numerical Mathematica implementation of these formulas. 
When the formulas of our paper and that of ref. |44| are both implemented with appropriate 
caching in C++, for the NMHV tree amplitudes for Vqqggggg and VqqQQggg, their evaluation 
time is similar |47|. 

We remark that in approaches based on BCFW recursion relations, the asymptotic number 
of terms in N fc MHV n amplitudes as n becomes large is quadratic in n for NMHV, quartic for 
NNMHV and worse for higher k. This is the reason we especially simplified the NMHV and 
NNMHV formulas presented in our paper, since we expect that they will be the most useful 
for practical applications, especially for small n. For k > 2 and large n there are at least two 
efficient numerical strategies making use of these formulae. First, one could use our formulae as 
initial conditions for a numerical implementation of the BCFW recursion relations, as described 
in section 3. Alternatively, one could use the Berends-Giele approach for k > 2, implemented 
using an efficient caching, in combination with our formulas for k < 2 |48|. 

We also included further functions that evaluate directly the simplified amplitudes of Ap- 
pendix A and B, They can be accessed via the following functions. 
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GGTnmhvgluon [n, a, b] 

is the simplified n-parton NMHV gluon amplitude with negative-helicity gluons at positions 
a, b and n. 

GGTnnmhvgluon [n, a, b, c] 

is the simplified n-parton NNMHV gluon amplitude with negative-helicity gluons at posi- 
tions a, b, c and n. 

GGTnmhv2f erm [n, c, a, b~\ 

is the simplified n-parton NMHV two-fermion amplitude with negative-helicity gluons at 
positions c, n and a fermion/anti-fermion at positions a and b. 

GGTnnmhv2f erm [n, c\, c 2 , a, b~\ 

is the simplified n-parton NNMHV two-fermion amplitude with negative-helicity gluons at 
positions c±, c 2 and n and a fermion/anti-fermion at position a and b. 

GGTnmhv4f erm[n, {{a 1: A t }, {a 2 , A 2 }}, {{b u Si}, {6 2 , 5 2 }}] 

is the simplified n-parton NMHV four-fermion amplitude with a negative-helicity gluon at 
position n, two gluinos of flavors A; t at positions and two anti-gluinos of flavors B; L at 
positions 6j. 

GGTnmhv4f ermS [n, {ai, a 2 }, {61, 6 2 }] 

is the simplified n-parton NMHV four-fermion amplitude with a negative-helicity gluon at 
position n and equally flavored gluinos/anti-gluinos at positions and 6j, respectively. 

GGTnnmhv4f erm[n, c, {{ai, AJ, {a 2 , A 2 }}, {{61, Si}, {b 2 , B 2 }}] 

is the simplified n-parton NNMHV four-fermion amplitude with two negative-helicity gluons 
at position c, n, two gluinos of flavors at positions and two anti-gluinos of flavors Bi 
at positions bi. 

GGTnnmhv4f ermS [n, c, {ai, a 2 }, {61, 6 2 }] 

is the simplified n-parton NNMHV four-fermion amplitude with negative-helicity gluons at 
positions c, n and equally flavored gluinos/anti-gluinos at positions a« and bi, respectively. 

GGTnmhv6f erm [n, B 3 , {{a u A ± }, {a 2 , A 2 }, {a 3 , A 3 }}, {{b u B ± }, {b 2 , B 2 }}] 

is the simplified n-parton NMHV six-fermion amplitude with three gluinos of flavors A4 at 

positions Oj and three anti-gluinos of flavors Bi at positions bi. Note that b 3 = n. 

GGTnmhv6f ermS [n, {ai, a 2 , a 3 }, {bi, b 2 }1 

is the simplified n-parton NMHV six-fermion amplitude with equally flavored gluinos/anti- 
gluinos at positions a« and bi, respectively. Note that 63 = n. 

GGTnnmhv6f erm[n, {{a u A x }, {a 2 , A 2 }, {a 3 , A 3 }}, {{b u B x }, {b 2 , B 2 }, {b 3 , B 3 }}] 
is the simplified n-parton NNMHV six-fermion amplitude with a negative helicity gluon at 
position n and three gluinos of flavors Ai at positions a« and three anti-gluinos of flavors 
Bi at positions 6j. 
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• GGTnnmhv6f ermS In, {a%, a>2, a 3 }, {b\, b<i, ^3}] 

is the simplified n-parton NNMHV six-fermion amplitude with a negative helicity gluon at 
position n and equally flavored gluinos/anti-gluinos at positions Oj and hi, respectively. 

The full list of functions available in GGT can be obtained by typing 

$GGTf unctions 

along with the documentation of each implemented function that can be accessed via the com- 
mand 

?GGTgluon 

for example. 
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